3.1.54 \(\int \sin ^4(a+b x) \sqrt {d \tan (a+b x)} \, dx\) [54]

3.1.54.1 Optimal result
3.1.54.2 Mathematica [A] (verified)
3.1.54.3 Rubi [A] (warning: unable to verify)
3.1.54.4 Maple [B] (warning: unable to verify)
3.1.54.5 Fricas [C] (verification not implemented)
3.1.54.6 Sympy [F]
3.1.54.7 Maxima [A] (verification not implemented)
3.1.54.8 Giac [A] (verification not implemented)
3.1.54.9 Mupad [F(-1)]

3.1.54.1 Optimal result

Integrand size = 21, antiderivative size = 257 \[ \int \sin ^4(a+b x) \sqrt {d \tan (a+b x)} \, dx=-\frac {21 \sqrt {d} \arctan \left (1-\frac {\sqrt {2} \sqrt {d \tan (a+b x)}}{\sqrt {d}}\right )}{32 \sqrt {2} b}+\frac {21 \sqrt {d} \arctan \left (1+\frac {\sqrt {2} \sqrt {d \tan (a+b x)}}{\sqrt {d}}\right )}{32 \sqrt {2} b}+\frac {21 \sqrt {d} \log \left (\sqrt {d}+\sqrt {d} \tan (a+b x)-\sqrt {2} \sqrt {d \tan (a+b x)}\right )}{64 \sqrt {2} b}-\frac {21 \sqrt {d} \log \left (\sqrt {d}+\sqrt {d} \tan (a+b x)+\sqrt {2} \sqrt {d \tan (a+b x)}\right )}{64 \sqrt {2} b}-\frac {7 \cos ^2(a+b x) (d \tan (a+b x))^{3/2}}{16 b d}-\frac {\cos ^4(a+b x) (d \tan (a+b x))^{7/2}}{4 b d^3} \]

output
-21/64*arctan(1-2^(1/2)*(d*tan(b*x+a))^(1/2)/d^(1/2))*d^(1/2)/b*2^(1/2)+21 
/64*arctan(1+2^(1/2)*(d*tan(b*x+a))^(1/2)/d^(1/2))*d^(1/2)/b*2^(1/2)+21/12 
8*ln(d^(1/2)-2^(1/2)*(d*tan(b*x+a))^(1/2)+d^(1/2)*tan(b*x+a))*d^(1/2)/b*2^ 
(1/2)-21/128*ln(d^(1/2)+2^(1/2)*(d*tan(b*x+a))^(1/2)+d^(1/2)*tan(b*x+a))*d 
^(1/2)/b*2^(1/2)-7/16*cos(b*x+a)^2*(d*tan(b*x+a))^(3/2)/b/d-1/4*cos(b*x+a) 
^4*(d*tan(b*x+a))^(7/2)/b/d^3
 
3.1.54.2 Mathematica [A] (verified)

Time = 0.49 (sec) , antiderivative size = 122, normalized size of antiderivative = 0.47 \[ \int \sin ^4(a+b x) \sqrt {d \tan (a+b x)} \, dx=-\frac {\left (21 \arcsin (\cos (a+b x)-\sin (a+b x)) \csc (a+b x) \sqrt {\sin (2 (a+b x))}+21 \csc (a+b x) \log \left (\cos (a+b x)+\sin (a+b x)+\sqrt {\sin (2 (a+b x))}\right ) \sqrt {\sin (2 (a+b x))}+18 \sin (2 (a+b x))-2 \sin (4 (a+b x))\right ) \sqrt {d \tan (a+b x)}}{64 b} \]

input
Integrate[Sin[a + b*x]^4*Sqrt[d*Tan[a + b*x]],x]
 
output
-1/64*((21*ArcSin[Cos[a + b*x] - Sin[a + b*x]]*Csc[a + b*x]*Sqrt[Sin[2*(a 
+ b*x)]] + 21*Csc[a + b*x]*Log[Cos[a + b*x] + Sin[a + b*x] + Sqrt[Sin[2*(a 
 + b*x)]]]*Sqrt[Sin[2*(a + b*x)]] + 18*Sin[2*(a + b*x)] - 2*Sin[4*(a + b*x 
)])*Sqrt[d*Tan[a + b*x]])/b
 
3.1.54.3 Rubi [A] (warning: unable to verify)

Time = 0.43 (sec) , antiderivative size = 251, normalized size of antiderivative = 0.98, number of steps used = 14, number of rules used = 13, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.619, Rules used = {3042, 3071, 252, 252, 266, 826, 1476, 1082, 217, 1479, 25, 27, 1103}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \sin ^4(a+b x) \sqrt {d \tan (a+b x)} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \sin (a+b x)^4 \sqrt {d \tan (a+b x)}dx\)

\(\Big \downarrow \) 3071

\(\displaystyle \frac {d \int \frac {(d \tan (a+b x))^{9/2}}{\left (\tan ^2(a+b x) d^2+d^2\right )^3}d(d \tan (a+b x))}{b}\)

\(\Big \downarrow \) 252

\(\displaystyle \frac {d \left (\frac {7}{8} \int \frac {(d \tan (a+b x))^{5/2}}{\left (\tan ^2(a+b x) d^2+d^2\right )^2}d(d \tan (a+b x))-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 252

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{4} \int \frac {\sqrt {d \tan (a+b x)}}{\tan ^2(a+b x) d^2+d^2}d(d \tan (a+b x))-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 266

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{2} \int \frac {d^2 \tan ^2(a+b x)}{d^4 \tan ^4(a+b x)+d^2}d\sqrt {d \tan (a+b x)}-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 826

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{2} \left (\frac {1}{2} \int \frac {d^2 \tan ^2(a+b x)+d}{d^4 \tan ^4(a+b x)+d^2}d\sqrt {d \tan (a+b x)}-\frac {1}{2} \int \frac {d-d^2 \tan ^2(a+b x)}{d^4 \tan ^4(a+b x)+d^2}d\sqrt {d \tan (a+b x)}\right )-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 1476

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{2} \left (\frac {1}{2} \left (\frac {1}{2} \int \frac {1}{d^2 \tan ^2(a+b x)-\sqrt {2} d^{3/2} \tan (a+b x)+d}d\sqrt {d \tan (a+b x)}+\frac {1}{2} \int \frac {1}{d^2 \tan ^2(a+b x)+\sqrt {2} d^{3/2} \tan (a+b x)+d}d\sqrt {d \tan (a+b x)}\right )-\frac {1}{2} \int \frac {d-d^2 \tan ^2(a+b x)}{d^4 \tan ^4(a+b x)+d^2}d\sqrt {d \tan (a+b x)}\right )-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 1082

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{2} \left (\frac {1}{2} \left (\frac {\int \frac {1}{-d^2 \tan ^2(a+b x)-1}d\left (1-\sqrt {2} \sqrt {d} \tan (a+b x)\right )}{\sqrt {2} \sqrt {d}}-\frac {\int \frac {1}{-d^2 \tan ^2(a+b x)-1}d\left (\sqrt {2} \sqrt {d} \tan (a+b x)+1\right )}{\sqrt {2} \sqrt {d}}\right )-\frac {1}{2} \int \frac {d-d^2 \tan ^2(a+b x)}{d^4 \tan ^4(a+b x)+d^2}d\sqrt {d \tan (a+b x)}\right )-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 217

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{2} \left (\frac {1}{2} \left (\frac {\arctan \left (\sqrt {2} \sqrt {d} \tan (a+b x)+1\right )}{\sqrt {2} \sqrt {d}}-\frac {\arctan \left (1-\sqrt {2} \sqrt {d} \tan (a+b x)\right )}{\sqrt {2} \sqrt {d}}\right )-\frac {1}{2} \int \frac {d-d^2 \tan ^2(a+b x)}{d^4 \tan ^4(a+b x)+d^2}d\sqrt {d \tan (a+b x)}\right )-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 1479

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{2} \left (\frac {1}{2} \left (\frac {\int -\frac {\sqrt {2} \sqrt {d}-2 \sqrt {d \tan (a+b x)}}{d^2 \tan ^2(a+b x)-\sqrt {2} d^{3/2} \tan (a+b x)+d}d\sqrt {d \tan (a+b x)}}{2 \sqrt {2} \sqrt {d}}+\frac {\int -\frac {\sqrt {2} \left (\sqrt {d}+\sqrt {2} \sqrt {d \tan (a+b x)}\right )}{d^2 \tan ^2(a+b x)+\sqrt {2} d^{3/2} \tan (a+b x)+d}d\sqrt {d \tan (a+b x)}}{2 \sqrt {2} \sqrt {d}}\right )+\frac {1}{2} \left (\frac {\arctan \left (\sqrt {2} \sqrt {d} \tan (a+b x)+1\right )}{\sqrt {2} \sqrt {d}}-\frac {\arctan \left (1-\sqrt {2} \sqrt {d} \tan (a+b x)\right )}{\sqrt {2} \sqrt {d}}\right )\right )-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{2} \left (\frac {1}{2} \left (-\frac {\int \frac {\sqrt {2} \sqrt {d}-2 \sqrt {d \tan (a+b x)}}{d^2 \tan ^2(a+b x)-\sqrt {2} d^{3/2} \tan (a+b x)+d}d\sqrt {d \tan (a+b x)}}{2 \sqrt {2} \sqrt {d}}-\frac {\int \frac {\sqrt {2} \left (\sqrt {d}+\sqrt {2} \sqrt {d \tan (a+b x)}\right )}{d^2 \tan ^2(a+b x)+\sqrt {2} d^{3/2} \tan (a+b x)+d}d\sqrt {d \tan (a+b x)}}{2 \sqrt {2} \sqrt {d}}\right )+\frac {1}{2} \left (\frac {\arctan \left (\sqrt {2} \sqrt {d} \tan (a+b x)+1\right )}{\sqrt {2} \sqrt {d}}-\frac {\arctan \left (1-\sqrt {2} \sqrt {d} \tan (a+b x)\right )}{\sqrt {2} \sqrt {d}}\right )\right )-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{2} \left (\frac {1}{2} \left (-\frac {\int \frac {\sqrt {2} \sqrt {d}-2 \sqrt {d \tan (a+b x)}}{d^2 \tan ^2(a+b x)-\sqrt {2} d^{3/2} \tan (a+b x)+d}d\sqrt {d \tan (a+b x)}}{2 \sqrt {2} \sqrt {d}}-\frac {\int \frac {\sqrt {d}+\sqrt {2} \sqrt {d \tan (a+b x)}}{d^2 \tan ^2(a+b x)+\sqrt {2} d^{3/2} \tan (a+b x)+d}d\sqrt {d \tan (a+b x)}}{2 \sqrt {d}}\right )+\frac {1}{2} \left (\frac {\arctan \left (\sqrt {2} \sqrt {d} \tan (a+b x)+1\right )}{\sqrt {2} \sqrt {d}}-\frac {\arctan \left (1-\sqrt {2} \sqrt {d} \tan (a+b x)\right )}{\sqrt {2} \sqrt {d}}\right )\right )-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

\(\Big \downarrow \) 1103

\(\displaystyle \frac {d \left (\frac {7}{8} \left (\frac {3}{2} \left (\frac {1}{2} \left (\frac {\arctan \left (\sqrt {2} \sqrt {d} \tan (a+b x)+1\right )}{\sqrt {2} \sqrt {d}}-\frac {\arctan \left (1-\sqrt {2} \sqrt {d} \tan (a+b x)\right )}{\sqrt {2} \sqrt {d}}\right )+\frac {1}{2} \left (\frac {\log \left (-\sqrt {2} d^{3/2} \tan (a+b x)+d^2 \tan ^2(a+b x)+d\right )}{2 \sqrt {2} \sqrt {d}}-\frac {\log \left (\sqrt {2} d^{3/2} \tan (a+b x)+d^2 \tan ^2(a+b x)+d\right )}{2 \sqrt {2} \sqrt {d}}\right )\right )-\frac {(d \tan (a+b x))^{3/2}}{2 \left (d^2 \tan ^2(a+b x)+d^2\right )}\right )-\frac {(d \tan (a+b x))^{7/2}}{4 \left (d^2 \tan ^2(a+b x)+d^2\right )^2}\right )}{b}\)

input
Int[Sin[a + b*x]^4*Sqrt[d*Tan[a + b*x]],x]
 
output
(d*(-1/4*(d*Tan[a + b*x])^(7/2)/(d^2 + d^2*Tan[a + b*x]^2)^2 + (7*((3*((-( 
ArcTan[1 - Sqrt[2]*Sqrt[d]*Tan[a + b*x]]/(Sqrt[2]*Sqrt[d])) + ArcTan[1 + S 
qrt[2]*Sqrt[d]*Tan[a + b*x]]/(Sqrt[2]*Sqrt[d]))/2 + (Log[d - Sqrt[2]*d^(3/ 
2)*Tan[a + b*x] + d^2*Tan[a + b*x]^2]/(2*Sqrt[2]*Sqrt[d]) - Log[d + Sqrt[2 
]*d^(3/2)*Tan[a + b*x] + d^2*Tan[a + b*x]^2]/(2*Sqrt[2]*Sqrt[d]))/2))/2 - 
(d*Tan[a + b*x])^(3/2)/(2*(d^2 + d^2*Tan[a + b*x]^2))))/8))/b
 

3.1.54.3.1 Defintions of rubi rules used

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 217
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(-(Rt[-a, 2]*Rt[-b, 2])^( 
-1))*ArcTan[Rt[-b, 2]*(x/Rt[-a, 2])], x] /; FreeQ[{a, b}, x] && PosQ[a/b] & 
& (LtQ[a, 0] || LtQ[b, 0])
 

rule 252
Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> Simp[c*(c*x 
)^(m - 1)*((a + b*x^2)^(p + 1)/(2*b*(p + 1))), x] - Simp[c^2*((m - 1)/(2*b* 
(p + 1)))   Int[(c*x)^(m - 2)*(a + b*x^2)^(p + 1), x], x] /; FreeQ[{a, b, c 
}, x] && LtQ[p, -1] && GtQ[m, 1] &&  !ILtQ[(m + 2*p + 3)/2, 0] && IntBinomi 
alQ[a, b, c, 2, m, p, x]
 

rule 266
Int[((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> With[{k = De 
nominator[m]}, Simp[k/c   Subst[Int[x^(k*(m + 1) - 1)*(a + b*(x^(2*k)/c^2)) 
^p, x], x, (c*x)^(1/k)], x]] /; FreeQ[{a, b, c, p}, x] && FractionQ[m] && I 
ntBinomialQ[a, b, c, 2, m, p, x]
 

rule 826
Int[(x_)^2/((a_) + (b_.)*(x_)^4), x_Symbol] :> With[{r = Numerator[Rt[a/b, 
2]], s = Denominator[Rt[a/b, 2]]}, Simp[1/(2*s)   Int[(r + s*x^2)/(a + b*x^ 
4), x], x] - Simp[1/(2*s)   Int[(r - s*x^2)/(a + b*x^4), x], x]] /; FreeQ[{ 
a, b}, x] && (GtQ[a/b, 0] || (PosQ[a/b] && AtomQ[SplitProduct[SumBaseQ, a]] 
 && AtomQ[SplitProduct[SumBaseQ, b]]))
 

rule 1082
Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> With[{q = 1 - 4*S 
implify[a*(c/b^2)]}, Simp[-2/b   Subst[Int[1/(q - x^2), x], x, 1 + 2*c*(x/b 
)], x] /; RationalQ[q] && (EqQ[q^2, 1] ||  !RationalQ[b^2 - 4*a*c])] /; Fre 
eQ[{a, b, c}, x]
 

rule 1103
Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> S 
imp[d*(Log[RemoveContent[a + b*x + c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, 
e}, x] && EqQ[2*c*d - b*e, 0]
 

rule 1476
Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[ 
2*(d/e), 2]}, Simp[e/(2*c)   Int[1/Simp[d/e + q*x + x^2, x], x], x] + Simp[ 
e/(2*c)   Int[1/Simp[d/e - q*x + x^2, x], x], x]] /; FreeQ[{a, c, d, e}, x] 
 && EqQ[c*d^2 - a*e^2, 0] && PosQ[d*e]
 

rule 1479
Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[ 
-2*(d/e), 2]}, Simp[e/(2*c*q)   Int[(q - 2*x)/Simp[d/e + q*x - x^2, x], x], 
 x] + Simp[e/(2*c*q)   Int[(q + 2*x)/Simp[d/e - q*x - x^2, x], x], x]] /; F 
reeQ[{a, c, d, e}, x] && EqQ[c*d^2 - a*e^2, 0] && NegQ[d*e]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3071
Int[sin[(e_.) + (f_.)*(x_)]^(m_)*((b_.)*tan[(e_.) + (f_.)*(x_)])^(n_.), x_S 
ymbol] :> With[{ff = FreeFactors[Tan[e + f*x], x]}, Simp[b*(ff/f)   Subst[I 
nt[(ff*x)^(m + n)/(b^2 + ff^2*x^2)^(m/2 + 1), x], x, b*(Tan[e + f*x]/ff)], 
x]] /; FreeQ[{b, e, f, n}, x] && IntegerQ[m/2]
 
3.1.54.4 Maple [B] (warning: unable to verify)

Leaf count of result is larger than twice the leaf count of optimal. \(618\) vs. \(2(197)=394\).

Time = 13.87 (sec) , antiderivative size = 619, normalized size of antiderivative = 2.41

method result size
default \(\frac {\left (16 \sqrt {2}\, \sqrt {-\frac {\cos \left (b x +a \right ) \sin \left (b x +a \right )}{\left (\cos \left (b x +a \right )+1\right )^{2}}}\, \left (\cos ^{3}\left (b x +a \right )\right ) \sin \left (b x +a \right )+16 \left (\cos ^{2}\left (b x +a \right )\right ) \sin \left (b x +a \right ) \sqrt {2}\, \sqrt {-\frac {\cos \left (b x +a \right ) \sin \left (b x +a \right )}{\left (\cos \left (b x +a \right )+1\right )^{2}}}-44 \cos \left (b x +a \right ) \sin \left (b x +a \right ) \sqrt {2}\, \sqrt {-\frac {\cos \left (b x +a \right ) \sin \left (b x +a \right )}{\left (\cos \left (b x +a \right )+1\right )^{2}}}-44 \sin \left (b x +a \right ) \sqrt {2}\, \sqrt {-\frac {\cos \left (b x +a \right ) \sin \left (b x +a \right )}{\left (\cos \left (b x +a \right )+1\right )^{2}}}+21 \ln \left (-\frac {\cot \left (b x +a \right ) \cos \left (b x +a \right )-2 \cot \left (b x +a \right )-2 \sin \left (b x +a \right ) \sqrt {-\left (\cot ^{3}\left (b x +a \right )\right )+3 \left (\cot ^{2}\left (b x +a \right )\right ) \csc \left (b x +a \right )-3 \cot \left (b x +a \right ) \left (\csc ^{2}\left (b x +a \right )\right )+\csc ^{3}\left (b x +a \right )+\cot \left (b x +a \right )-\csc \left (b x +a \right )}-2 \cos \left (b x +a \right )-\sin \left (b x +a \right )+\csc \left (b x +a \right )+2}{-1+\cos \left (b x +a \right )}\right )-21 \ln \left (-\frac {\cot \left (b x +a \right ) \cos \left (b x +a \right )-2 \cot \left (b x +a \right )+2 \sin \left (b x +a \right ) \sqrt {-\left (\cot ^{3}\left (b x +a \right )\right )+3 \left (\cot ^{2}\left (b x +a \right )\right ) \csc \left (b x +a \right )-3 \cot \left (b x +a \right ) \left (\csc ^{2}\left (b x +a \right )\right )+\csc ^{3}\left (b x +a \right )+\cot \left (b x +a \right )-\csc \left (b x +a \right )}-2 \cos \left (b x +a \right )-\sin \left (b x +a \right )+\csc \left (b x +a \right )+2}{-1+\cos \left (b x +a \right )}\right )-42 \arctan \left (\frac {-\sin \left (b x +a \right ) \sqrt {2}\, \sqrt {-\frac {\cos \left (b x +a \right ) \sin \left (b x +a \right )}{\left (\cos \left (b x +a \right )+1\right )^{2}}}+\cos \left (b x +a \right )-1}{-1+\cos \left (b x +a \right )}\right )+42 \arctan \left (\frac {\sin \left (b x +a \right ) \sqrt {2}\, \sqrt {-\frac {\cos \left (b x +a \right ) \sin \left (b x +a \right )}{\left (\cos \left (b x +a \right )+1\right )^{2}}}+\cos \left (b x +a \right )-1}{-1+\cos \left (b x +a \right )}\right )\right ) \sqrt {d \tan \left (b x +a \right )}\, \cos \left (b x +a \right ) \sqrt {2}}{128 b \left (\cos \left (b x +a \right )+1\right ) \sqrt {-\frac {\cos \left (b x +a \right ) \sin \left (b x +a \right )}{\left (\cos \left (b x +a \right )+1\right )^{2}}}}\) \(619\)

input
int(sin(b*x+a)^4*(d*tan(b*x+a))^(1/2),x,method=_RETURNVERBOSE)
 
output
1/128/b*(16*2^(1/2)*(-cos(b*x+a)*sin(b*x+a)/(cos(b*x+a)+1)^2)^(1/2)*cos(b* 
x+a)^3*sin(b*x+a)+16*cos(b*x+a)^2*sin(b*x+a)*2^(1/2)*(-cos(b*x+a)*sin(b*x+ 
a)/(cos(b*x+a)+1)^2)^(1/2)-44*cos(b*x+a)*sin(b*x+a)*2^(1/2)*(-cos(b*x+a)*s 
in(b*x+a)/(cos(b*x+a)+1)^2)^(1/2)-44*sin(b*x+a)*2^(1/2)*(-cos(b*x+a)*sin(b 
*x+a)/(cos(b*x+a)+1)^2)^(1/2)+21*ln(-(cot(b*x+a)*cos(b*x+a)-2*cot(b*x+a)-2 
*sin(b*x+a)*(-cot(b*x+a)^3+3*cot(b*x+a)^2*csc(b*x+a)-3*cot(b*x+a)*csc(b*x+ 
a)^2+csc(b*x+a)^3+cot(b*x+a)-csc(b*x+a))^(1/2)-2*cos(b*x+a)-sin(b*x+a)+csc 
(b*x+a)+2)/(-1+cos(b*x+a)))-21*ln(-(cot(b*x+a)*cos(b*x+a)-2*cot(b*x+a)+2*s 
in(b*x+a)*(-cot(b*x+a)^3+3*cot(b*x+a)^2*csc(b*x+a)-3*cot(b*x+a)*csc(b*x+a) 
^2+csc(b*x+a)^3+cot(b*x+a)-csc(b*x+a))^(1/2)-2*cos(b*x+a)-sin(b*x+a)+csc(b 
*x+a)+2)/(-1+cos(b*x+a)))-42*arctan((-sin(b*x+a)*2^(1/2)*(-cos(b*x+a)*sin( 
b*x+a)/(cos(b*x+a)+1)^2)^(1/2)+cos(b*x+a)-1)/(-1+cos(b*x+a)))+42*arctan((s 
in(b*x+a)*2^(1/2)*(-cos(b*x+a)*sin(b*x+a)/(cos(b*x+a)+1)^2)^(1/2)+cos(b*x+ 
a)-1)/(-1+cos(b*x+a))))*(d*tan(b*x+a))^(1/2)*cos(b*x+a)/(cos(b*x+a)+1)/(-c 
os(b*x+a)*sin(b*x+a)/(cos(b*x+a)+1)^2)^(1/2)*2^(1/2)
 
3.1.54.5 Fricas [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.40 (sec) , antiderivative size = 947, normalized size of antiderivative = 3.68 \[ \int \sin ^4(a+b x) \sqrt {d \tan (a+b x)} \, dx=\text {Too large to display} \]

input
integrate(sin(b*x+a)^4*(d*tan(b*x+a))^(1/2),x, algorithm="fricas")
 
output
1/256*(16*(4*cos(b*x + a)^3 - 11*cos(b*x + a))*sqrt(d*sin(b*x + a)/cos(b*x 
 + a))*sin(b*x + a) + 21*b*(-d^2/b^4)^(1/4)*log(9261/2*d^2*cos(b*x + a)*si 
n(b*x + a) + 9261/2*(b^3*(-d^2/b^4)^(3/4)*cos(b*x + a)^2 - b*d*(-d^2/b^4)^ 
(1/4)*cos(b*x + a)*sin(b*x + a))*sqrt(d*sin(b*x + a)/cos(b*x + a)) - 9261/ 
4*(2*b^2*d*cos(b*x + a)^2 - b^2*d)*sqrt(-d^2/b^4)) - 21*b*(-d^2/b^4)^(1/4) 
*log(9261/2*d^2*cos(b*x + a)*sin(b*x + a) - 9261/2*(b^3*(-d^2/b^4)^(3/4)*c 
os(b*x + a)^2 - b*d*(-d^2/b^4)^(1/4)*cos(b*x + a)*sin(b*x + a))*sqrt(d*sin 
(b*x + a)/cos(b*x + a)) - 9261/4*(2*b^2*d*cos(b*x + a)^2 - b^2*d)*sqrt(-d^ 
2/b^4)) + 21*I*b*(-d^2/b^4)^(1/4)*log(9261/2*d^2*cos(b*x + a)*sin(b*x + a) 
 - 9261/2*(I*b^3*(-d^2/b^4)^(3/4)*cos(b*x + a)^2 + I*b*d*(-d^2/b^4)^(1/4)* 
cos(b*x + a)*sin(b*x + a))*sqrt(d*sin(b*x + a)/cos(b*x + a)) + 9261/4*(2*b 
^2*d*cos(b*x + a)^2 - b^2*d)*sqrt(-d^2/b^4)) - 21*I*b*(-d^2/b^4)^(1/4)*log 
(9261/2*d^2*cos(b*x + a)*sin(b*x + a) - 9261/2*(-I*b^3*(-d^2/b^4)^(3/4)*co 
s(b*x + a)^2 - I*b*d*(-d^2/b^4)^(1/4)*cos(b*x + a)*sin(b*x + a))*sqrt(d*si 
n(b*x + a)/cos(b*x + a)) + 9261/4*(2*b^2*d*cos(b*x + a)^2 - b^2*d)*sqrt(-d 
^2/b^4)) + 21*b*(-d^2/b^4)^(1/4)*log(9261*d^2 + 18522*(b^3*(-d^2/b^4)^(3/4 
)*cos(b*x + a)*sin(b*x + a) - b*d*(-d^2/b^4)^(1/4)*cos(b*x + a)^2)*sqrt(d* 
sin(b*x + a)/cos(b*x + a))) - 21*b*(-d^2/b^4)^(1/4)*log(9261*d^2 - 18522*( 
b^3*(-d^2/b^4)^(3/4)*cos(b*x + a)*sin(b*x + a) - b*d*(-d^2/b^4)^(1/4)*cos( 
b*x + a)^2)*sqrt(d*sin(b*x + a)/cos(b*x + a))) + 21*I*b*(-d^2/b^4)^(1/4...
 
3.1.54.6 Sympy [F]

\[ \int \sin ^4(a+b x) \sqrt {d \tan (a+b x)} \, dx=\int \sqrt {d \tan {\left (a + b x \right )}} \sin ^{4}{\left (a + b x \right )}\, dx \]

input
integrate(sin(b*x+a)**4*(d*tan(b*x+a))**(1/2),x)
 
output
Integral(sqrt(d*tan(a + b*x))*sin(a + b*x)**4, x)
 
3.1.54.7 Maxima [A] (verification not implemented)

Time = 0.41 (sec) , antiderivative size = 225, normalized size of antiderivative = 0.88 \[ \int \sin ^4(a+b x) \sqrt {d \tan (a+b x)} \, dx=\frac {21 \, d^{6} {\left (\frac {2 \, \sqrt {2} \arctan \left (\frac {\sqrt {2} {\left (\sqrt {2} \sqrt {d} + 2 \, \sqrt {d \tan \left (b x + a\right )}\right )}}{2 \, \sqrt {d}}\right )}{\sqrt {d}} + \frac {2 \, \sqrt {2} \arctan \left (-\frac {\sqrt {2} {\left (\sqrt {2} \sqrt {d} - 2 \, \sqrt {d \tan \left (b x + a\right )}\right )}}{2 \, \sqrt {d}}\right )}{\sqrt {d}} - \frac {\sqrt {2} \log \left (d \tan \left (b x + a\right ) + \sqrt {2} \sqrt {d \tan \left (b x + a\right )} \sqrt {d} + d\right )}{\sqrt {d}} + \frac {\sqrt {2} \log \left (d \tan \left (b x + a\right ) - \sqrt {2} \sqrt {d \tan \left (b x + a\right )} \sqrt {d} + d\right )}{\sqrt {d}}\right )} - \frac {8 \, {\left (11 \, \left (d \tan \left (b x + a\right )\right )^{\frac {7}{2}} d^{6} + 7 \, \left (d \tan \left (b x + a\right )\right )^{\frac {3}{2}} d^{8}\right )}}{d^{4} \tan \left (b x + a\right )^{4} + 2 \, d^{4} \tan \left (b x + a\right )^{2} + d^{4}}}{128 \, b d^{5}} \]

input
integrate(sin(b*x+a)^4*(d*tan(b*x+a))^(1/2),x, algorithm="maxima")
 
output
1/128*(21*d^6*(2*sqrt(2)*arctan(1/2*sqrt(2)*(sqrt(2)*sqrt(d) + 2*sqrt(d*ta 
n(b*x + a)))/sqrt(d))/sqrt(d) + 2*sqrt(2)*arctan(-1/2*sqrt(2)*(sqrt(2)*sqr 
t(d) - 2*sqrt(d*tan(b*x + a)))/sqrt(d))/sqrt(d) - sqrt(2)*log(d*tan(b*x + 
a) + sqrt(2)*sqrt(d*tan(b*x + a))*sqrt(d) + d)/sqrt(d) + sqrt(2)*log(d*tan 
(b*x + a) - sqrt(2)*sqrt(d*tan(b*x + a))*sqrt(d) + d)/sqrt(d)) - 8*(11*(d* 
tan(b*x + a))^(7/2)*d^6 + 7*(d*tan(b*x + a))^(3/2)*d^8)/(d^4*tan(b*x + a)^ 
4 + 2*d^4*tan(b*x + a)^2 + d^4))/(b*d^5)
 
3.1.54.8 Giac [A] (verification not implemented)

Time = 0.31 (sec) , antiderivative size = 245, normalized size of antiderivative = 0.95 \[ \int \sin ^4(a+b x) \sqrt {d \tan (a+b x)} \, dx=\frac {\frac {42 \, \sqrt {2} {\left | d \right |}^{\frac {3}{2}} \arctan \left (\frac {\sqrt {2} {\left (\sqrt {2} \sqrt {{\left | d \right |}} + 2 \, \sqrt {d \tan \left (b x + a\right )}\right )}}{2 \, \sqrt {{\left | d \right |}}}\right )}{b} + \frac {42 \, \sqrt {2} {\left | d \right |}^{\frac {3}{2}} \arctan \left (-\frac {\sqrt {2} {\left (\sqrt {2} \sqrt {{\left | d \right |}} - 2 \, \sqrt {d \tan \left (b x + a\right )}\right )}}{2 \, \sqrt {{\left | d \right |}}}\right )}{b} - \frac {21 \, \sqrt {2} {\left | d \right |}^{\frac {3}{2}} \log \left (d \tan \left (b x + a\right ) + \sqrt {2} \sqrt {d \tan \left (b x + a\right )} \sqrt {{\left | d \right |}} + {\left | d \right |}\right )}{b} + \frac {21 \, \sqrt {2} {\left | d \right |}^{\frac {3}{2}} \log \left (d \tan \left (b x + a\right ) - \sqrt {2} \sqrt {d \tan \left (b x + a\right )} \sqrt {{\left | d \right |}} + {\left | d \right |}\right )}{b} - \frac {8 \, {\left (11 \, \sqrt {d \tan \left (b x + a\right )} d^{5} \tan \left (b x + a\right )^{3} + 7 \, \sqrt {d \tan \left (b x + a\right )} d^{5} \tan \left (b x + a\right )\right )}}{{\left (d^{2} \tan \left (b x + a\right )^{2} + d^{2}\right )}^{2} b}}{128 \, d} \]

input
integrate(sin(b*x+a)^4*(d*tan(b*x+a))^(1/2),x, algorithm="giac")
 
output
1/128*(42*sqrt(2)*abs(d)^(3/2)*arctan(1/2*sqrt(2)*(sqrt(2)*sqrt(abs(d)) + 
2*sqrt(d*tan(b*x + a)))/sqrt(abs(d)))/b + 42*sqrt(2)*abs(d)^(3/2)*arctan(- 
1/2*sqrt(2)*(sqrt(2)*sqrt(abs(d)) - 2*sqrt(d*tan(b*x + a)))/sqrt(abs(d)))/ 
b - 21*sqrt(2)*abs(d)^(3/2)*log(d*tan(b*x + a) + sqrt(2)*sqrt(d*tan(b*x + 
a))*sqrt(abs(d)) + abs(d))/b + 21*sqrt(2)*abs(d)^(3/2)*log(d*tan(b*x + a) 
- sqrt(2)*sqrt(d*tan(b*x + a))*sqrt(abs(d)) + abs(d))/b - 8*(11*sqrt(d*tan 
(b*x + a))*d^5*tan(b*x + a)^3 + 7*sqrt(d*tan(b*x + a))*d^5*tan(b*x + a))/( 
(d^2*tan(b*x + a)^2 + d^2)^2*b))/d
 
3.1.54.9 Mupad [F(-1)]

Timed out. \[ \int \sin ^4(a+b x) \sqrt {d \tan (a+b x)} \, dx=\int {\sin \left (a+b\,x\right )}^4\,\sqrt {d\,\mathrm {tan}\left (a+b\,x\right )} \,d x \]

input
int(sin(a + b*x)^4*(d*tan(a + b*x))^(1/2),x)
 
output
int(sin(a + b*x)^4*(d*tan(a + b*x))^(1/2), x)